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ABSTRACT 

We present the source catalog and the properties of the B-, R-, and /—band 
images obtained to support the AKARI North Ecliptic Pole Wide (NEP-Wide) 
survey. The NEP-Wide is an AKARI infrared imaging survey of the north ecliptic 
pole covering a 5.8 deg 2 area over 2.5 - 6 fim wavelengths. The optical imaging 
data were obtained at the Maidanak Observatory in Uzbekistan using the Seoul 
National University 4k x 4k Camera on the 1.5m telescope. These images cover 
4.9 deg 2 where no deep optical imaging data are available. Our B— , R— , and 
/—band data reach the depths of ~23.4, ~23.1, and ~22.3 mag (AB) at 5a, 
respectively. The source catalog contains 96,460 objects in the R— band, and 
the astrometric accuracy is about 0.15" at la in each RA and Dec direction. 
These photometric data will be useful for many studies including identification 
of optical counterparts of the infrared sources detected by AKARI, analysis of 
their spectral energy distributions from optical through infrared, and the selection 
of interesting objects to understand the obscured galaxy evolution. 

Subject headings: catalogs — galaxies: evolution — galaxies: photometry - 
surveys 
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INTRODUCTION 



Recently, the observational study of the cosmic star formation history has been at the 
center of our efforts to understand the evolution and the formation of galaxies. This is be- 
cause the observational data play a key role in constraining physical processes governing the 
evolutionary history of galaxies on which outcomes of the theoretical predictions are strongly 



galaxies (e.g., 


Goto et al. 


2010; 


Le Borgne et al. 


2008; 


Shim et al. 


2007; 


Im 


2005 


), however, there 



the 



dust extinction of the light coming from star formation. 



Since the discovery of ultra-luminous infrared galaxies (ULIRGs) and luminous infrared 
galaxies (LIRGs) by the InfraRed Astronomical Satellite (IRAS), it has been recognized that 
a significant portion of the light from the star formation is missing in the optical/ultraviolet 
but emitted in infrared (IR). Therefore, statistical studies of infrared luminous galaxies to 
reveal the hidden cosmic star formation activities were major themes of the more recent 
infrared space telescope missions such as the Infrared Space Observatory (ISO), the Spitzer 
Space Telescope, and the AKARI Space Telescope. Such studies found t hat the star for- 
mation at higher redshifts (z ~ 1) are more obscured than at z ~ (e.g., iLe Floc'h et al. 



20051 ; iBell et al.l 120071 ; iGoto et al.l 120101 ) underlining the importance of infrared-based stud- 
ies of galaxies. However, due to the lack of efficient coverage between 10 to 20 ^m in the 
Spitzer imaging, comparison of the same rest-frame mid-infrared (MIR) quantities over a 
wide redshift range is not possible to further examine the usefulness of the MIR star forma- 
tion indicators at various redshifts and the correlation between the star formation and the 
relative strengths of MIR spectroscopic features made by molecules such as Poly Aromatic 
Hydrocarbon (PAH) and silicates. The additional MIR coverage at 10 t o 20 nm is also 



advantageous in studying the MIR-excess emission in early-type galaxies (IKo et al. 



2007 



Shim et al.l 120101). and the identification of Active Galactic Nuclei (AGN; e.g., iLee et al. 



2009 



Huang et al.l 120091 ). 



The AKARI satellite is an infrared space telescope launched in February 2006. By 
offering a contiguous, wide-field imaging capability in the wavelength range between 2.5 
and 26/im, AKARI complements the capability of the Spitzer space telescope, and provides 
an unprecedented view of the universe, especially at 11 and 15/mi. Because of the Sun- 
synchronous orbit of AKARI, deep observations can be done only at the North and the 
South ecliptic poles. Therefore, the North Ecliptic Pole (NEP) survey is one of the important 
extragalactic surveys of AKARI. The NEP survey has two major programs, NEP- Wide and 
NEP-Deep surveys. The NEP- Wide field covers 5.8 deg 2 centered at NEP, whereas the NEP- 
Deep field covers 0.38 deg 2 with deeper exposures than the NEP- Wide. Due to the wide field, 
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the contiguous MIR w avelength coverage, and the low extinction value of E(B — V) = 0.047 
at the center of NEP ( ISchlegel et al.lll998l ). the NEP-Wide field is unique for the studies of 
the extragalactic objects like infrared luminou s galaxies and AGNs. Fo r further informatio n 
abo ut NEP surveys, the r eaders can refer to iMatsuhara et al.l (120061 ) , IWada et al.l ( 120071 ) , 
and lLee et all (120071 Eooih . 



Optical imaging data can enhance the value of the unique AKARI NEP survey. Since 
AKARI has a low spatial resolution (the full width at half maximum (FWHM) of the point 
spread function (PSF) is about 6" at 15/xm), the optical imaging with a higher resolution 
can help improving the astrometric accuracy and help deblending confused objects. Also, 
the optical bands are important to construct multi-wavelength spectral energy distributions 
(SEDs) of infrared sources such as galaxies and AGNs to understand their properties. There 
already exist deep optical data covering 2 deg 2 of the central region of the NEP-Wide field, 
obtained by the Ca nada-France-Hawaii Telescope (CFHT) with u*,g',r',i', and z' filters 
( iHwang et al.l 120071 1 . However, there is no deep optical data for the remaining NEP-Wide 
area. The Digitized Sky Survey images are available, but neither their depths nor spatial 
resolutions are sufficient (the depth is about 20 mag at R and the seeing is about 3" with 
pixel scale of 1.0" at the NEP area) for identifying optical counterparts of many infrared 
sources. Hence, we observed this remaining NEP-Wide area in optical B, R, and I filters. 
The BRI filters were chosen because we needed a multiple color set at least three colors to 
construct the optical part of the SED of galaxy; one blue part, one red part, and a middle 
part between red and infrared bands. Also, those Bessell filters are wider than the Sloan 
Digital Sky Survey (SDSS) filters, enabling us to achieve a higher signal to noise ratio (S/N) 
per unit exposure time than the SDSS filters. 

Throughou t the paper, we use a cosmology with Q M = 0.3, Q\ = 0.7, and H = 70 km 



s 1 Mpc 1 (e.g.. llm et al.lll997l ). We also use the AB magnitude system. 



2. OBSERVATION 

The AKARI NEP-Wide field is centered at a = 18 h 00 m 00 s , 5 = +66°36'00" covering 
a circular area of 5.8 deg 2 . Our optical BRI imaging survey covers 4.9 deg 2 outside the 
central 2 deg 2 area of the AKARI NEP-Wide. Figure [1] shows the field coverage of different 
components of the NEP-Wide survey. The thick, solid circle indicates the AKARI NEP- 
Wide (5.8 deg 2 ) coverage. The dashed squares are the CFHT optical survey area. The small 
squares show the Maidanak optical imaging area. It consists of 70 subfields, each of them 
having an 18' x 18' field of view. These Maidanak fields are grouped into 4 different sections, 
East (8 fields), West (7 fields), North (23 fields), and South (32 fields), and we name each 
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subfield with E01-E08, W01-W08, N01-N23, and S01-S32. The names of the subfields and 
their central positions are listed in Table [TJ 

The images were obtained in Bessell B, R, and I filters from 2007 June 12 to August 
5 at the Maidanak Observatory in Uzbekistan, using the Seoul Nation al University 4k x 4k 



Camera (SNUCAM) on the 1.5 m Richey- Chretien, AZT-22 telescope dim et all 120101 ). The 
CCD camera has a 4096 x 4096 format and a pixel scale of 0.266" at f/7.75. We used the 
four amplifier mode for the CCD image readout. 

Twilight flats were taken at the start and the end of each night when possible. Since the 
/—band images are found to have complicated fringe patterns, all of the /—band images were 
dithered by ~20" in random directions. These dithered images were used to construct the 
mast er fringe i mage . For photometric calibration, we observed photometric standard fields 



from lLandoltl (119921 ) at two different altitudes at least. In most subfields, standard star 
images were obtained on the same night, but in a few cases no corresponding standard star 
observations were carried out due to bad weather. In such cases, the photometric calibration 
was done using a different way, as described in section 13.31 

Table [2] summarizes our observations. For each filter, typical total exposure times are 
21 min in B, 30 min in R, and 20 min in / filter. The images for a given subfield and a given 
filter were obtained mostly on a single night so that images are in homogeneous quality at 
least for a single filter of a given subfield. There are a few exceptions to this rule (3 subfields 
for B— band and 13 subfields for /—band). This happened when we needed to take images 
of a given subfield and a given filter over two nights due to observational constraints such as 
the lack of the observing time (e.g., the night was ending before finishing all the images of 
the subfield), or a rapidly deteriorating weather or seeing. In some cases, the total exposure 
times and the number of images of each subfield were adjusted to the weather conditions. 
For example, when there were thin clouds in the sky, the exposure time or the number of 
images were increased to compensate for the loss of the S/N. The frame times for individual 
exposures were also adjusted. Shorter exposure times (2 minutes) were used when there 
were no guide stars, or the weather condition was varying too rapidly. The B— and R— band 
images were mostly taken during the same nights under dark/grey sky to provide images 
with similar seeing values, but the /—band images were taken mostly during the bright time 
since the depth of /—band images is the least influenced by the presence of the moon among 
the BRI filters. The seeing condition for each subfield was determined as the median value 
of FWHM of non-saturated point sources with the magnitudes between 13 and 15 mag. In 
general, the seeing FWHM values of the stacked R— band and /—band images are stable at 
~ 1.2" and ~ 1-1", respectively, but there were more variations in the seeing conditions in 
B— band and the B— band images have the worst seeing FWHM values in general among the 
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BRI filters. The maximum and the minimum seeing values among the BRI filter images 
of a given subfield differ by a factor of 1.5 in median, but in extreme cases, the difference 
could be as large as a factor of ~2 due to deteriorated seeing conditions in the B— band. The 
5a detection limits of each subfield were determined from the sky background sigma, within 
an aperture diameter of 3 times FWHM value. The detection limits vary from subfield to 
subfield, but the median 5a detection limit in R— band is 23.1 mag. Table [3] lists the seeings 
and 5a depths for each subfield. Figure [2] shows the seeing and the 5a detection limit of the 
three filters for each subfield. 



3. DATA REDUCTION 

3.1. Preprocessing 

We used IRAF0 noao . imred. ccdred package to process raw images. The data were 
processed with usual reduction procedures such as bias subtractions using a combined bias 
image and flat-field correction using twilight flat of each night. Dark correction was unnec- 
essary since there is virtually no dark current in the images taken with our CCD camera. 

Because the bias subtraction and the flat fielding could not remove the amplifier bias 
level variations perfectly due to small temporal variations in the bias level, some of the pre- 
processed images show slightly different bias levels in different quadrants, each representing 
a section of image read by a single amplifier. To remove this pattern, we determined the 
background values of the 4 quadrants, and subtracted the differences between the back- 
ground levels of the 4 quadrant images from each quadrant. Then the background levels are 
normalized to the background level of a quadrant with the lowest background value. 

We removed the fringe pattern in J— band following the steps below. First, we created 
a master fringe image combining all /—band images (716 images) in median. Then, before 
subtracting the master fringe frame, we re-scaled the master fringe frame by an appropriate 
factor. This scaling factor is taken as the ratio of the median background value of each 
science image to the master fringe frame background value. Finally, we subtracted the re- 
scaled master fringe frame from each /—band image. We checked the difference in image 
quality before and after the fringe removal. Before the fringe pattern removal, the large scale 
fringe pattern varies ±3% of the background level which is comparable to the sky sigma value 



1 IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Associa- 
tion of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 
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per pixel. After the fringe subtraction, the large scale variation is reduced to < 0.5% of the 
background or more than 6 times less than the sky sigma value, showing that the fringe 
pattern is removed successfully well below the sky sigma values. 

After the fringe correction and the registration, all images were combined in median. 
Cosmic rays are removed during the stacking of multiple images with a sigma clipping adopt- 
ing +3-sigma for rejecting outlier pixels. 



3.2. Astrometry 

To derive astrometry solutions, we used a software named SCAMP which finds astrom- 
etry solution by cross-matching automatically an input reference catalog containing a set of 
accurate RA and Dec values of sources against a list of object position s obtained from an 



input FITS image ( lBertinll2006l ). In our case, we used USNO-B1 catalog ( iMonet et al.ll2003l ) 
as the reference catalog for the astrometry calibration. Astrometry solutions were derived 
for all the stacked, subfield images in each filter separately. The degree of the polynomial for 
distortion correction is adopted to be 4, since it gave the best solution among other options. 
SCAMP can refine the astrometry solution by comparing positions of objects in overlapping 
images, but the procedure does not apply to our case. 

To check the accuracy of the astrometry solution, we compared our catalog with the 
USNO-B1 catalog stars that had the positional errors of < 0.2" and the apparent magnitudes 
of 14 < Bl mag < 19. Figure [3]- (a) shows the distribution of the differences in position 
between the USNO-B1 stars and the corresponding objects in our catalog in RA and Dec 
directions. Figures [3]- (b) and (c) show the histograms of the astrometric differences between 
our catalog and USNO-B1 catalog of RA and Dec, respectively. To derive the astrometry 
error, we fit the histogram of the astrometric differences with Gaussian functions. The dotted 
lines indicate the best-fit Gaussian function to the distributions. After the fitting, we find 
that the la astrometric errors are 0.16" for RA and 0.15" for Dec. 



3.3. Photometric Calibration 

We carried out the photometry calibration using three methods depending on t he avail- 



abilit y of the standard star data. The first method uses the Landolt standard stars (JLandolt 



19921 ). This procedure applies to 25 fields in B— band, 25 fields in R— band, and 18 fields 
in /—band, where the standard star data were taken at two different altitudes at the least. 
The adopted standard calibration formula is rrii ns — M st d = ko + k\X, where rrii ns is an in- 
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strumental magnitude ( rrii ns = 25 — 2.5 log (DN/ sec) ), M s ^ is a standard magnitude, ko is 
a zero point, k\ is an airmass coefficient, and X is an airmass. The instrumental magnitudes 
were measured within an aperture of a 12" diameter on the standard stars. We checked with 
growth curves that this aperture is large enough to contain virtually all photons from the star 
irrespective of the seeing of our data. We chose the standard stars that have similar colors 
for deriving the coefficients of the above formula. To confirm the effect of not including the 
color terms, we checked the difference between the listed magnitudes of standard stars with 
various colors (those not included in the derivation of the photometric zero-points) and the 
magnitudes of these standard stars derived with our photometry solutions. We find that the 
rms difference is less than 0.05 mag. That is, the magnitudes derived with our photometry 
solutions are independent of the color terms well within the error 0.05 mag. This is consis- 
tent with the expected variation in the photometric zero points due to neglecting the color 
term, base d on an observat ional campaign of photometry standards using the SNUCAM at 



Maidanak (ILim et al.l 120091 ) . Therefore we decided to neglect the color term of the formula. 
We found ko and k\ values for each night and each filter, and derived the standard magni- 
tudes of the objects using these values. The ko has the values in the range of 1.67 ± 0.13 at 
B— band, 1.70 ± 0.11 at R— band, and 2.07 ± 0.08 at J— band. Also, the k\ have the values 
of 0.33 ± 0.05 at 5-band, 0.10 ± 0.03 at i?-band, and 0.04 ± 0.01 at /-band. T hese values 



are c onsistent with values derived from a campaign observation of standard stars (ILim et al. 



20091 ). The zero-point errors were estimated from the rms magnitude differences between 
the known magnitudes of the standard stars and the magnitudes of the same standard stars 
derived from our zero-point solution. The zero-point error is about 0.032 mag using this 
method. 

The second photometric calibration method is adopted for the data taken during the 
nights which have only one altitude data for standard stars (34, 40, and 27 fields at B-, R-, 
and /—band, respectively). In such cases, we used other night's k\ value since the night-to- 
night variation in k\ is not significant. This produces an additional zero-point error of 0.028 
mag, at most. Consequently, the total error associated with the zero-point determination of 
the second photometric calibration method is less than 0.040 mag. 

The third method of the photometric calibration is applied to the subfields without 
associated standard star observations. For such fields (11, 5, and 25 fields in B-, R-, and 
/—band, respectively), we used photometry of sources in the overlapping area of a neighbor 
subfield. The corresponding sources that have magnitudes between 12 and 16 mag were used 
to measure the photometric zero-point. For this method, the zero-point error was derived 
from the standard deviation of the magnitude offset between the sources in the field where 
we want to determine the photometric zero-point and the matched sources in the neighbor 
fields. In this case, the zero-point error is about 0.045 mag with respect to the photometry 
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of the neighbor field. Therefore, the total error in the zero-point for the third method is 
about 0.032-0.045 mag, depending on the method used for the photometry calibration of 
the neighbor field. 

We checked the consistency of our photometric calibration using an overlap area between 
two neighboring fields (those that were not used for the photometric calibration of the other 
field). In the case of the first method, the mean photometry differences are — 0.006± 0.028 
mag in B (30 fields for overlaps), 0.002 ± 0.033 mag in R (23 fields for overlaps), and 0.001 
± 0.017 mag in I magnitudes (17 fields for overlaps). In the case of the second method, 
the mean photometry differences are about 0.002 ± 0.030 mag in B (37 fields for overlaps), 
— 0.006± 0.040 mag in R (44 fields for overlaps), and — 0.022± 0.035 mag in / magnitudes 
(20 fields for overlaps). Note that the consistency of our photometric calibration between 
different nights are well below 0.05 mag for all filters which is expected from the zero-point 
errors. 

The Z p for each subfield are listed in Table |3j For the nights which have the ko and 
the ki values, the Z p means 25 — (k + kiX), i.e. the zero point for converting DN per 
exposure time to magnitude, including the airmass correction. Therefore we can determine 
the apparent magnitude M, through M = Z p — 2.5 log (DN/sec). The Z p can be used for 
the purpose of checking the depth of each subfield. Figure 0] shows several depths of each 
subfield for the three filters. For each subfield, the depth are checked using the 5a detection 
limit, the 50% completeness (see section FOj) . the 99% reliability (see section I4T5]) . or the Z p , 
as well as the total exposure time. 



CATALOGS 



4.1. Object Detection and Photometry 



SExtractor ( IBertin &: Arnoutslll996l ) was used for the detection and the photometry of 
objects in the images. We stacked B-, R-, and /—band images together, and used the BRI 
stacked image as our object detection image. After many trials of SExtractor parameter sets 
and examining the detected and the non-detected sources by eye, we chose a set of 1.2cr for 
DETECT_THRESH, 5 connected pixels for DETECT_MINAREA, DEBLEND_NTHRESH 
of 64, DEBLEND_MINCONT of 0.005, 200 for BACK_SIZE, and 3 for BACK_FILTERSIZE 
as the optimal set of detection parameters. Lowering the detection limit produced too many 
spurious detections and increasing the detection limit tend to miss obvious sources. The 
adopted detection limit corresponds to the S/N of abou t 4.5 which is usually an optimal 
S/N limit for detection of objects (e.g., IShim et al.l 120061 ). 
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To cross- identify and perform the photometry of objects in each band, we used the 
ASSOC parameter in SExtractor with a matching radius of 0.658", 3a value of the astrometric 
error. The above process may miss faint objects which are brighter in a particular filter than 
the others. To supplement the above source detection, we also ran SExtractor on images 
in each filter separately. This process reveals faint sources detected in a single filter image, 
but not detected in the BRI stacked image. These sources detected in a single filter image 
were also matched with detections in other single filter images with the matching radius 
of 0.658". Adding these faint sources, we detected 63,333 for B, 96,460 for R and 70,492 
sources for / filter over the reliability of 99% (see section I4~5l) . Among these objects, 16% 



in B, 33% in R, and 15% in / are detected only in the single filter images. We derived 
aperture magnitudes using an apert ure of diam eter of 3 times FWHM and auto magnitudes 



with Kron-like elliptical apertures (lKronlll980l ) which are taken as the total magnitudes 



The aperture correction value is 0.10 ± 0.02 mag for all the filters and all the subfields. 

The magnitude errors are computed by combining the standard SExtractor estimates 
and the zero-point errors we estimated. The standard SExtractor error estimates are based 
on the Poisson statistics, and could underestimate the noise values when th e sky noises in 



adjoi ning pixels are correlated due to the st acking of the d i there d frames (IGawiser et al. 



20061 ). Using a similar procedure described in IGawiser et al.l (120061 ). we calculated how the 



error within a square aperture of N x N pixels, ctjv, changes as a function of N. In the 
case where each pixel noise is uncorrelated with each other, ctjv oc N. When the noise is 
correlated, the exponent of N is greater than 1. In our data, we find oc A^ 1 ' 14 . For a 
magnitude measured within an aperture with 3.0" diameter (3 times a typical FWHM of our 
image), there are about 110 pixels (or N ~ 10). Therefore, the photometric errors of faint 
objects (R > 22 mag) could be underestimated by as much as about a factor of 1.4. For the 
brighter objects (R < 20 mag), the photometric errors are not affected by this effect, since 
their errors are dominated by the object noise. 

For the same object that exist on two or more fields in the overlapping areas, we chose 
an object that had the smallest FLAGS value from SExtractor (the FLAGS has a nonzero 
value if the object is blended, saturated, and/or truncated at an edge of a image) and the 
smallest photometric error. Among the same objects, we only included the one with a better 
S/N or image quality in the final catalog. The examples of such cases are objects on the 
edge of an image or objects elongated due to astigmatism near the edge of the CCD in one 
image, but do not suffer such a problem in the other. At the bright end, the photometry is 
free from saturation and non-linearity above R ~ 12 mag. We cut off the faint sources over 
the magnitude of false detection rate 1% (see section 143]) . 



We found spurious detections due to the reflection halos caused by the bright objects 
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and due to faint, uncorrected fringe patterns in R— band images. In case of the reflection 
halos around bright stars, we divided the bright stars into three classes, and set a radius 
to define the region where spurious detections lie within it (120" at 8 < B mag < 10 , 80" 
at 10 < B mag < 12 , 40" at 12 < B mag < 14). We flag detections inside these radii as 
'near_bright_obj', since they are likely to be spurious detections or their photometry are not 
reliable. We also found many spurious detections in the outer regions of stacked R— band 
images, thus we flag objects that are not matched with B— or /—band catalogs and exist 
within 30" from the edge as 'near_edge'. We found 5,186 (5%) of 'near_bright_obj' at R— band 
and 14,833 (15%)of 'near_edge' sources, and checked those spurious sources by eyes. These 
flags are included in the BRI merged catalog. Also we flag non-stellar sources as 'galaxy' 
(see, section I4T2]) . The flags are marked at the last column of the catalog (Table Hj). 



4.2. Star-Galaxy Separation 

To distinguish stellar sources from other extragalactic sources, we used a B — H vs. 
H — N2 color-color diagram and the SExtractor stellarity index. Here, the stellarity indices 
measure the likelihood of a source being a point source or an extended source, with the 
value of 1 for a perfect point source, and for a diffuse, extended source. The stellarities 
are measured from the single-filter stacked images. However, when we determine whether a 
source is point-like or extended, we adopted the stellarity value measured in a stacked image 
in the filter where the photometric error of the source is the smallest. Note that the stellarity 
values are nearly identical among different filter images when an object is sufficiently bright 
(R < 19 mag). We also note that the stellarity values are determined from the R— band 
images for most of the objects, since the R— band images generally have the highest S/N 
among the three filters. 

The if— band magnitudes come from the catalog made from the J— and if— band 
images of the NEP-Wide field which were took with FLAMINGOS on the 2.1m telescope 
at Kitt Peak National Observatory (Jeon et a!., in preparation), and the N2 magnitudes 



come from the AKARI observation, of which the effective wavelength is 2.43/im (ILee et al. 



20071 ). We matched the objects detected in H— and N2— band with the objects detected in 
B-, R-, and /—band using a matching radius of 1.5" (for H— band) and 3" (for N2— band), 
respectively. 



Figure [5] shows the B — H vs. H — N2 diagram. iBarkhouse and Halll (120011 ) used B — J 
vs. J — K colors to distinguish quasars from stars. Instead of J — K color, we used a similar 
selection method using H — N2. On the left panel of Figure the dots and the contours 
are made from the object detected in B— band with the magnitude error less than 0.1 and 



- 11 - 



matched with H and N2 detections. We can notice that there is a stellar locus at H — N2 
~ —0.5 mag with < B — H < 5. Because most of stars have similar slopes at wavelength 
greater than lfim on their SEDs, stars have H — N2 ~ —0.5 mag. We checked the stellar 
locus using stars (the crosses of the right panel of Figure [5]) in NEP-Wide field which had 
the stellarity greater than 0.8 and the magnitude of 10 < H < 13. Most of the stars are 
located at the stellar locus with a few exceptions. 

We also checked the location of ga laxies using redsh ift tracks of galaxies. The dashed 



line is for a star-forming galaxy, M51 (jSilva et al.l Il998[ ) and the dotted line is for a pas 



sively evolving galaxy (see figure caption for the mode l used for this). We considered the 



intergalactic medium (IGM) attenuation (lMadaulll996f ) for all redshift tracks. Most of the 



galaxy track lie in the area of H — N2 greater than —0.2 mag, but galaxies with z < 0.1 lie 
at H — N2 < —0.2 mag. The PSF of galaxy at low redshift (z < 0.1) is in general larger than 
that of star, so such sources can be identified as extended sources in the image. However, 
if the low redshift galaxies are faint and small, they may get misclassified as stars with the 
color cut. 

Since the magnitude limit of the H and N2— band are different from those in the optical, 
this method can be applied only to a relatively bright magnitude limit. Also, in order to 
identify bright, extended sources (namely, low redshift galaxies) in the stellar region on the 
color-color plot, we apply a stellarity cut where stellarity index has a value between 1 (point 
sources) and (extended sources). 

In summary, we define stellar sources as the following: 

Objects with H — N2 < —0.2 mag and the stellarity index > 0.2 or 
Objects with H — N2 > —0.2 mag and the stellarity index > 0.8 or 
if there is no H — N2 color available, stellarity index > 0.8 

The first condition weeds out extended objects at low redshift, while the second condition 
allows us to select clear, point sources in the region of galaxies on the color-color magnitude. 
The third criterion applies to a small number of objects near the detection limit of optical 
images, since H and N2 images show detections of 98% and 77% of objects brighter than 
the magnitude cut of the catalog at 50% completeness limit. 



4.3. Catalog Format 



The B-, R-, and /—band merged catalog is presented in Table HI The catalog contains 
objects that are deemed to be reliable above 99% confidence where the confidence values 
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are derived independently for each subfield in each filter (see section 14. 5p . Our catalog 
contains each object's identification number (ID), RA, Dec, the total magnitude and its 
error in each band, the aperture magnitude within a circular aperture with a diameter of 3 
times FWHM and its error in each band, the stellarity value fr om SExtractor, the galactic 
extinction value from the extinction map of ISchlegel et al.l (119981 ) for each band, and the flags 
described in section 14.11 Note that the listed magnitudes are not corrected for the galactic 
extinction. The conversion formulas of Vega to AB magnitudes are S(Vega) = _B(AB) + 0.09, 
i?(Vega) = R(AB) - 0.22 and /(Vega) = /(AB) - 0.45. T he offset values between V ega and 
AB magnitude are computed fro m the Vega s pectrum (IBohlin &: Gillilandl 120041 ) and B, 
R, and / filter response functions (|Besselllll990l ) coupled with the CCD quantum efficiency 
curve. The magnitude error is the square root of the quadratic summation of the photometric 
measurement error and the zero-point error from the photometric calibration. We use a 
dummy value of 99.00 for non-detection. 

Figure [6] shows the stellarity distribution (the left three panels) and the photometric 
error distribution (the right three panels) for all sources detected in B-, R-, and /— band 
images. The white lines on the right panels show the mode values for each magnitude bins, 
with the 25% and 75% quartiles indicated with error bars. The spread of the data points 
in the magnitude versus error plots are caused mainly by varying depths of the data from 
different subfields. On average, the photometric error is less than 0.1 mag at B < 22.5 mag, 
R < 22 mag, and / < 21 mag. 



4.4. Completeness 

Since the source detection near the detection limit is dominated by the detections per- 
formed on the single filter image (~ 90% of the faintest sources of R— band are detected 
only in the single filter image) rather than the BRI stacked images, we computed the com- 
pleteness using the B-, R-, and /—band single filter images of each subfield, to examine 
the completeness of the source detection. To estimate the completeness, we created artifi- 
cial objects for each magnitude bin from 12 to 23.5 mag and placed them randomly in the 
background regions of the single filter images. The artificial objects are made of 100 point 
sources and 100 extended sources where the elliptical galaxies account for 40% of the ex- 
tended sources. These artificial objects are created using the noao . artdata package of the 
IRAF. After performing detection and photometry of these artificial sources with the same 
parameter setting as those applied to the actual images, we calculated the completeness, 
which is defined as the ratio of the number of detected artificial sources to the number of 
input objects for each magnitude bin. Figure [7] shows the completeness of B—,R—, and /— 
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band data. The gray solid lines in Figure [7] represent the completeness of all the subfields 
calculated with this method. We overploted the deepest (S05 in B, S05 in R and N23 in 
I) and the shallowest (N22 in B, N18 in R, and S15 in J) fields in the 5a depths with the 
black solid lines. The galaxy number counts (see section 15. ip after the corrections based 
on these completeness result match the number counts from deeper surveys, meaning that 
our completeness estimates are reasonable. Figure 0] shows the 50% completeness as the 
asterisks for each subfield, as well as the 5a depth, the 99% reliability (see section H~5j) . the 
Z p , and the total exposure time. The magnitudes of completeness 50% for each subfield are 
presented in Table [3j We can see that the magnitudes at 50% completeness are brighter 
than the ones of 5a detection limits by a few tenths of magnitude. 

4.5. Reliability 

We examined the reliability of our detections using false detection rates for each filter 
and for each subfield. Here, we define false detection rate as the ratio of the number of 
objects detected in a negative image to the number of objects detected in the original image 
at a given magnitude bin. The negative image means the original image multiplied by — 1. 
By this definition, an image with noises only would have the false detection rate of 1, while 
an image containing objects but with no noise would have the false detection rate of 0. 
When estimating the false detection rates, we excluded areas where objects are marked with 
'near_edge', since the inclusion of such area could disproportionately overestimate the false 
detection rate in the other areas which occupy most of the image. 

Figure El shows the false detection rates of B-, R-, and /—band data. The magnitude 
of the false detection rates 1% (or the reliability 99%) of all the subfields are indicated in 
Table [31 Also, Figure 0] shows the 99% reliability as diamond symbols for each subfield, as 
well as the 5a depth, 50% completeness, the Z p , and the total exposure time. Figure 0] show 
that the magnitude limit of reliability 99% is greater than the 50 % completeness limit for 
B— and R— band, except ~20 fields. The 99% reliability limit in /—band is consistently 
brighter than the 50% completeness limit by 0.3 mag. The systematic trend arises from 
false detections of the fringe patterns left in /—band image even after the fringe subtraction. 
Field-to-field variations of false detection rates and the completeness limits are much larger 
for .R-band than the other filters. It is due to the large difference of the image depths of each 
subfield in R— band, where some fields (from N02 to N23) have relatively small depth with 
the 5a limit < 22.5. It is because of the moon phase since these images were taken during 
full moon phase (29 July - 1 August). 

As a separate exercises, we examined the magnitude distribution of objects marked as 
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'near_edge' or 'near_bright_obj'. Note that most of these objects are likely to be spurious 
detections. 80% of these sources are found to be fainter than the magnitude limit of the 99% 
reliability. The 99% reliability cut of the catalog removes most of the spurious detections, 
but we note that objects marked with 'near_edge' or 'near_bright_obj' should be treated with 
caution. 



PROPERTIES OF THE DATA 



5.1. Galaxy Number Counts 



To obtain the depth and homogeneity of this survey, we show the galaxy number count 
as a function of magnitude. The total area in which we used the galaxy number count of 
bright galaxies is about 3.8 deg 2 where H— and N2— band data are also available. We take 
the Poisson error as the error estimate of number counts. 

We calculated the number of galaxies per square degree at each magnitude bin. Since 
the depth of each subfield is different from another, the combined number counts of all the 
NEP survey fields were computed following the steps below. First, the number count was 
computed for each subfield up to the completeness limit of 50% as determined in section 
14.41 Then, the completeness correction was applied to the number count of each subfield. 
Finally, we combined the number counts from all the subfields to derive the final number 
count result. With this method, the bright end of the number counts, usually susceptible 
to field-to-field variations, include the data from all the 70 fields, while the deepest number 
counts include the data from the deepest subfields only. We also calculated the number 
count of stars using the stellar objects of magnitude brighter than the completeness limit of 
~ 50%. 

Figure M displays the galaxy and the stellar number counts in B— , R— , and J— band. 
The lower solid line with crosses represents raw number count and the upper solid line 
with crosses is the number count after applying correction for the completeness. The 
solid line with asterisks shows the number count of star s. Al so we compare our num- 
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The galaxies and stars number count data are provided in Table [5j 



5.2. Color-Magnitude Diagram 



Figure [10] displays the color-magnitude diagrams of R vs. B — R (left) and its histograms 
(right). The upper panel shows the diagram for galaxies with R magnitude error less than 
0.1, and the lower panel shows for the stars as the same conditions with galaxies in the upper 
panel. We checked positions of galaxies on the color-magnitude diagram traced by the red- 



shifts . The dotted line represents a redshift track of a star-forming galaxy, M51 (jSilva et al. 



19981 ). and the dashed line is for a passively evolving galaxy. The SEP of a passively evolv- 
ing galaxy is for a 5 Gyr old stellar population, created with the iBruzual fc Charlotl (120031 ) 
model, assuming a spontaneous burst of single stellar population with metallicity of Z=0.02 
and the Salpeter IMF . We included the attenuation of galaxy light due to the IGM when 
computing the colors (IMadaul 119961 ) . Two different lines represent galaxies with character- 
istic absolute magnitudes: for an early type galaxy, the characteristic absolute magnitude, 
M*(r*), is taken as —20.75 + 5 l og h, and for a late type galaxy, —20.30 + 5 log h is adopted for 
M*(r*) (INakamura et al.l 120031 ). To convert the r* — ban d magnitude to the R— band magni- 
tude, we used the magnitude transformation relation of I Jester et al.l (120051 ). and g — r ~ 0.8 
as the mean color of galaxies. The k corrections are calculated from the templates of two 
galaxies in a standard manner. The B — R colors of elliptical and spiral galaxies are 1-1.5 
mag at z=0.1. For 0.2 < z < 0.6, the elliptical galaxies have the B — R values of 1.5-3 mag 
and the spiral galaxies have the B — R values of 1-1.5 mag. The color-magnitude plot shows 
that our B and R— band images detect galaxies at z = through z = 0.6, and galaxies 
brighter than M*(r*) to higher redshifts {z ~ 1). On the other hand, we find bimodality in 
the color-magnitude diagrams of stellar objects (bottom plot, Figure [TU]). The bimodality 
originates from G dwarfs that dominate the color at B — R ~ 0.5 and M Giant stars at B — R 
~ 2. Figure HU displays the color-magnitude diagrams in / vs. R — I (left) and its histograms 
(right) with the sources of the I magnitude error less than 0.1. The symbols and the tracks 
are the same as in Figure [JUJ The R — I color of elliptical and spiral galaxies is 0.3-0.5 at z 
< 0.2 and 0.5-0.9 at z ~ 0.5. Like Figure [T0~| the distribution of objects in the I versus R — I 
plot confirm that J— band images have the depth similar to the R— band images in terms of 
detecting galaxies at different redshifts. We can notice that the model colors of B — R and 
R — I are consistent with the distributions of galaxies. We expect that the combination of 
the optical data and the AKARI infrared data to provide a unique view on the evolution 
of galaxies, such as the morpho logical transform ation process through studying MID excess 
emission in early-type galaxies ( IKo et al.l 120091 ). 
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Discussion 



Here, we investigate how our data compare with other similar optical datasets, and 
explore the potential of selecting high redshift galaxies and quasars using the data. 

First, many extragalactic surveys have been performed with various depths and area 
coverage. The 5.8 deg 2 coverage o f the NEP- Wide sur vey is wider than that of the Cos- 
mic Evol ution Survey (CO SMOS, IScoville et al.l 120071 ). or other deeper surveys such as 



DEEP2 (Davis et al. 



20031 ) and the Great Observatories Origins Deep Survey (GOODS, 



Dickinson et al.ll2003l ). On the other hand, the area coverage is smaller than that of the 
Spitzer Wide- Area Infrared Extragalactic Survey (SWIRE) survey (ILonsdale et al.l 120031 ) by 
a factor of 10 with a similar depth in infrared as well as optical ancillary datasets. In terms 
of the dep th and area cover age, our survey is most similar to the Spitzer First Look Sur- 
vey (FLS, iFadda et al.l I2004J ) . Our infrared depth is a bit shallower than those of the FLS, 
while we provide a contiguous wavelength coverage including 11 to 15 micron. Consequently, 
our optical imaging depths are designed to match the infrared survey depths, thus about a 
magnitude shallower than the R— band depth of the FLS ancillary data. 

Next, we discuss the usefulness of our datasets for selecting high redshift objects such 
as Lyman Break Galaxies (LBGs) and quasars. Figure [T5] shows color-color diagrams in 
B — R vs. R — I. We plotted sources with R and / magnitude errors less than 0.07 mag as 
gray points which correspond to objects with R < 21.0-23.0 mag, and / < 20.4-21.8 mag. 
The dashed line with diamonds represents a redshift track for a star- forming galaxy (M51, 
Silva et al.lll998l ) and the dotted line with squares is a reds hift track of a median composite 
spectrum of SDSS quasars from IVanden Berk et al.l (120011 ). The att enuatio n of th e source 
light below Lyman a due to IGM is calculated using the model of iMadaul (119961 ). Using 
the redshift tracks as a guide, we define a selection bo x in the color-color diagram where 



objects at z > 3.5 can be found (B— band dr opouts; e.g.. ISteidel et al.lll999l ; iGiavalisco et al. 



2004 ; iKang fc Imll2009l ; iGlikman et al.ll2010l ). The selection box is defined with the following 
equations: 



B-R > 2 

B-R > i(R-I) + 0.l 



(la) 
(lb) 



Although LBGs could be found in this selection box in principle, the expected number 
of the B— dropout g alaxies in our data is very low. Extrapolating the surface number density 
of z ~ 4 LBGs from ISteidel et al.l (119991 ) assuming an exponential decrease of the brightest 
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objects as in the Schechter function, we expect about one or less than one B— dropout galaxies 
at / < 21.5 mag. On the other hand, the number density could be higher by a factor of 
10 or so, if the bright end of the luminosity function flattens due to AG N contribution as 
seen in the lumin osity functions of [/—band dropout objects at z ~ 3 (e.g.. lHunt et al.ll2004j ; 



Shim et al.ll2007l ). In our data, we identify ~800 B— dropout candidates in the selection box, 



but we believe that most of these are interlopers such as stars and low redshift galaxies due 
to the large mismatch with the expected number. 



The expected number of quasars over the 4.9 deg 2 is < 10 at I < 20 mag ([Richards et al. 



2006 ). or 10 and 20 at R < 21.5 mag and 22.5 mag, respectively ( Glikman et al. 201ol ). After 



inspecting images of the .B-dropout candidates visually and examining their mult i- wavelength 
SEDs from the optical through the MIR ll/xm-band of AKARI, we identify at least several 
objects that have the SED shapes consistent with those of high redshift quasars. Judging 
from the shapes of the SEDs, most of the interlopers are low redshift galaxies or late-type 
stars. The selection of quasar candidates at z > 3.5 is ongoing, and we expect that the 
identification of z ~ 4 quasars at these magnitude limits will be valuab le for constraining 
the faint-end slope of quasars at high redshift (e.g., iGlikman et al.l 120101 ) . 



7. SUMMARY 

We have shown the characteristics of B—,R—, and J— band data at Maidanak Obser- 
vatory from the follow-up imaging observation of the AKARI NEP-Wide field. Using the 
SNUCAM on the 1.5m telescope, we covered 4.9 deg 2 with the 5er depth of ~23.4, ~23.1, 
and ~22.3 magnitudes (AB) at B—,R—, and /—band, respectively. We detected 63,333 
sources in B— band, 96,460 objects in R— band, and 70,492 in /—band. These data are now 
being used to identifying optical counterparts for accurate astrometry of infrared sources 
and deblending of confused sources. Also, our data provide a multi-wavelength coverage 
of the NEP-Wide field, enabling detailed SED analysis. Through these SEDs, we can get 
key physical properties such as stellar mass and photometric redshifts, and select interesting 
objects such as high redshift quasars. Spectroscopic follow-up surveys are being conducted 
on infrared luminous objects using the information gathered from this imaging survey. 
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Fig. 1. — Field coverage of the NEP-Wide optical survey. The center of the image corre- 
sponds to the NEP at a = 18 h 00 m 00 s , S = +66°36'00". Our BRI imaging data (4.9 deg 2 ) 
are indicated with small squares, while the thick, solid circle and the dashed squares indicate 
AKARI (5.8 deg 2 ) and CFHT (2 deg 2 ) coverage, respectively. The BRI images consist of 
70 fields. The names of the subfields and their central positions are listed in Table HJ 
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Fig. 2. — The upper and the lower panels show the seeing FWHMs and the 5a detection 
limits for each subfield, respectively (the diamonds for B, the triangles for R, and the squares 
for /—band). 
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Fig. 3. — Distribution (a) of the astrometry differences between our catalog and USNO-B1 
catalog, and the histograms of the differences in RA (b) and Dec (c) directions. In (b) 
and (c), the dotted lines indicate the best-fit Gaussian functions to the histograms of the 
distribution. From the Gaussian function, the la astrometric errors are 0.16" for RA and 
0.15" for Dec. 
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Fig. 4. — The figure shows several measures of the image quality and the depth of each 
subfield, such as the 5a depth (the cross symbols), the 50% completeness (the asterisks; see 
section 144]) . the 99% reliability (the diamonds; see section 1431) . the photometric zero point, 
Z p (the triangles), and the total exposure time (the squares). 
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Fig. 5. — B — H vs. H — N2 color-color diagrams and the density contours of sources 
detected in B, H, and N2. The points and the contours in both figures are made with 
objects that have B— band magnitude errors less than 0.1 mag. The dashed line is a redshift 
track of a star- forming galaxy (M5 1 ) and the dotted line is for a passively evolving galaxy 
made with the iBruzual fc Charlotl (120031 ) model of a passively evolving 5 Gyr-old galaxy 
with spontaneous burst, metallicity of Z=0.02, and the Salpeter IMF. We use the vertical 
line as a boundary that distinguishes stellar sources against galaxies and AGNs. In the right 
panel, we overlay stars in the NEP-Wide field (crosses) to indicate the effectiveness of the 
color-cut in separating stars from galaxies. 
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I magniiude [AB] I magniiude [AB] 

Fig. 6. — SExtractor stellarity distribution and the photometric error distribution for all 
sources detected in B— , Ft—, and I— band images, respectively. The stellarity index has a 
value between 1 (point sources) and (extended sources). In the right panel, we calculated 
the mode values of photometric error for each magnitude bin and let the 25% and 75% 
quartiles of the distributions as the error for each point. For B < 22.5 mag, R < 22 mag, 
and I < 21 mag, the photometric error is less than 0.1 mag. 
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Fig. 7. — Completeness for B— , R-, and J— band. The gray solid lines are the completeness 
of all subfields from the simulation method, especially for the deepest fields (S05 for B, S05 
for R, and N23 for 7;the upper black solid lines) and the shallowest fields (N22 for B, N18 
for R, and S15 for 7;the lower black solid lines) in the 5a depths. The horizontal dashed 
line is for 50% completeness cut. The magnitudes of completeness 50% for each subfield are 
shown in Table El 
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Fig. 8. — False detection rate (the ratio of the number of objects detected in the negative 
image to the number of objects detected in the original image for each magnitude bin) of all 
subfields in B-, R-, and /—band. The black solid lines are for the deepest or the shallowest 
fields as same as Figured The magnitudes of false detection rate 1% for each field are shown 
in Table |3j 
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Fig. 9. — Galaxy and stellar number counts of R-, and J— band (two solid lines with 
crosses). The error estimates of number counts are the Poisson errors. The lower solid line 
is for galaxy number count without completeness correction, and the upper solid line is from 
the galaxy number count corrected using the completeness. We compare our number counts 
with the number counts in the literature (other symbols). The solid lines with asterisks are 
for the number count of stars. The galaxy and stellar number count data are provided in 
Table E 
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Fig. 10. — The color-magnitude diagram of R vs. B — R (left) and histograms showing the 
distribution B — R colors (right). The upper panel shows the galaxies with R magnitude 
error less than 0.1. The dotted line is a redshift track for a star-forming galaxy (M51) and 
the dashed line is for a passively evolving galaxy. The redshift track is for galaxies with the 
characteristic absolute magnitudes (see the text). The lower panel shows the stars of the 
same conditions in the upper panel. 
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Fig. 11.— I vs. R — I color-magnitude diagrams (left) and the distribution of R — I 
colors (right) with the sources of / magnitude error less than 0.1. The upper panel is for 
the galaxies and the lower panel is for the stars. The symbols are the same as those in the 
Figure [10J 
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Fig. 12. — The color-color diagrams of B — R vs. R — I and the redshift tracks of a star- 
forming galaxy and a quasar. The gray points are sources with R and I magnitude errors 
less than 0.07. The dashed line with diamond points represents a track for M51, and the 
dotted line with squar es is for a redshift t rack made from a median composite spectrum 
of SDSS quasars from IVanden Berk et al.l (120011 ). The black boxed area in the upper-left 
corner defines the selection region for high redshift objects at z > 3.5. The asterisks show 
the sources in this area and the asterisks with arrows indicate objects with 2cr lower limits 
forB — R colors, when the object was not detected in B— band. 
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Table 1. Names of each subfield in Figure [T] and their central positions. 
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Table 2. Observation summary. 



Filter 


Total exposure time [min] 


Exposure time per image [min] 


Seeing ["] 


5cr depth [AB mag] 


50% Completeness [AB mag] 


B 


20-45 (21) 


2-15 


0.8-2.1 (1.4) 


22.56-24.34 (23.40) 


22.05-23.86 (22.98) 


R 


18-40 (30) 


2-10 


0.8-1.7 (1.2) 


22.07-24.14 (23.09) 


21.75-23.57 (22.60) 


I 


16-64 (20) 


2 


0.8-1.5 (1.1) 


21.55-22.90 (22.33) 


21.16-22.37 (21.85) 



Note. — 



The values in the parentheses are the median values. 



Table 3. Photometric zero-point, seeing, 5<r detection limit, magnitude of the completeness 50%, and magnitude of 

the reliability 99% for each subfield. 
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23.34 


24.12 


23.56 





.92 


23.92 


23.57 


23.86 


23.03 


1.50 


21.61 


21 


.28 


21.25 


S09 


23.00 


1.35 


23.93 


23.44 


24.18 


23.55 


1 


01 


23.91 


23.37 


23.72 


23.14 


1.06 


22.37 


22 


.09 


21.97 


S10 


23.00 


1.24 


24.09 


23.75 


23.96 


23.55 


1 


.08 


23.81 


23.23 


23.60 


23.03 


1.08 


21.81 


21 


.57 


21.55 



CO 



Table 3 — Continued 



B — band R— band J — band 



Field Z p a Seeing Depth b comp50%° rclia99% d Z p a Seeing Dcpth b comp50%° relia99% d Z p a Seeing Depth b comp50% c relia99% d 

[mag] ["] [mag] [mag] [mag] [mag] ["] [mag] [mag] [mag] [mag] ["] [mag] [mag] [mag] 



Sll 


23 


.00 


1 


.39 


24 


08 


23 


.71 


23 


78 


23 


48 





.99 


23 


.63 


23 


18 


23 


.30 


22 


73 


1 


.07 


21 


.58 


21 


28 


20.86 


S12 


23 


.00 


1 


.36 


24 


.00 


23 


.36 


24 


03 


23 


.48 





.92 


23 


80 


23 


.27 


23 


.70 


22. 


68 


1 


.10 


22 


.06 


21 


78 


20.73 


S13 


23 


.00 


1 


.37 


24 


04 


23 


.48 


23 


.50 


23 


.48 





.98 


23 


.77 


23 


04 


23 


.61 


22 


78 


1 


.09 


21 


.89 


21 


35 


19.30 


S14 


23 


00 


1 


.41 


23 


97 


23 


.62 


23 


.96 


23 


.45 


1 


.08 


23 


.67 


23 


.25 


23 


.17 


23. 


29 


1 


.04 


22 


.51 


22 


.08 


21.58 


S15 


22 


.98 


1 


.31 


23 


.81 


23 


.37 


23 


.97 


23 


.45 





.90 


23 


.73 


23 


.17 


23 


.50 


22. 


18 


1 


.18 


21 


.55 


21 


.21 


20.28 


S16 


22 


.97 


1 


.36 


23 


.91 


23 


.44 


23 


.61 


23 


.45 





.84 


23 


.82 


23 


.37 


23 


.50 


23. 


28 


1 


.09 


22 


.60 


22 


.24 


21.77 


S17 


22 


.98 


1 


.29 


23 


79 


23 


.41 


23 


.68 


23 


45 





.85 


23 


60 


23 


14 


23 


.45 


23. 


27 


1 


.10 


22 


.51 


21 


.86 


21.00 


S18 


22 


.98 


1 


.53 


23 


.63 


23 


.19 


23 


.98 


23 


45 





.92 


23 


.75 


23 


.23 


23 


.11 


23 


21 


1 


.11 


22 


.40 


22 


00 


21.51 


S19 


22 


.97 


1 


.59 


23 


.57 


23 


.11 


23 


.44 


23 


.44 


1 


.02 


23 


45 


23 


06 


23 


.35 


23. 


21 


1 


01 


22 


.40 


21 


98 


21.64 


S20 


22 


.97 


1 


.45 


23 


.71 


23 


.30 


23 


.55 


23 


.38 





.95 


23 


60 


23 


01 


23 


24 


23. 


16 


1 


32 


21 


.87 


21 


.36 


21.40 


S21 


22 


.96 


1 


.68 


23 


45 


22 


.96 


23 


.35 


23 


.37 





.93 


23 


.35 


22 


93 


23 


.15 


23. 


24 


1 


.39 


21 


.83 


21 


.32 


21.51 


S22 


23 


.05 


1 


.49 


23 


72 


23 


.26 


23 


.72 


23 


.37 





.91 


23 


.41 


22 


94 


22 


.91 


23. 


33 


1 


34 


21 


.93 


21 


51 


21.82 


S23 


23 


04 


1 


.14 


23 


89 


23 


.34 


23 


87 


23 


.38 


1 


.09 


23 


.50 


22 


92 


23 


51 


23. 


04 


1 


32 


21 


.74 


21 


.29 


21.30 


S24 


22 


91 


1 


.15 


23 


66 


23 


.26 


23 


98 


23 


.38 


1 


13 


23 


55 


23 


.22 


23 


.18 


23. 


35 


1 


.05 


22 


.29 


21 


75 


21.55 


S25 


22 


89 


1 


.06 


23 


72 


23 


31 


23 


.66 


23 


.37 


1 


.26 


23 


46 


23 


.03 


23 


.08 


23. 


36 


1 


.07 


22 


.12 


21 


77 


21.63 


S26 


22 


62 


1 


.58 


23 


29 


22 


.87 


23 


.01 


23 


38 


1 


17 


23 


.24 


22 


82 


22 


.74 


23. 


36 


1 


.07 


22 


.24 


21 


.94 


21.23 


S27 


22 


.83 


1 


.00 


23 


63 


23 


.23 


23 


.36 


23 


.47 


1 


.25 


23 


.62 


22 


.99 


23 


.40 


23. 


36 


1 


13 


22 


.23 


21 


.78 


21.04 


S28 


22 


.82 





.96 


23 


64 


23 


.28 


23 


.50 


23 


.33 


1 


.46 


23 


.21 


22 


59 


23 


.19 


23. 


36 


1 


.08 


22 


.25 


21 


95 


21.39 


S29 


22 


98 


1 


.97 


22 


.92 


22 


41 


23 


.57 


23 


.37 


1 


.29 


23 


19 


22 


.83 


23 


.08 


23. 


36 


1 


.11 


22 


.22 


21 


92 


21.26 


S30 


22 


.52 


1 


.88 


23 


02 


22 


.42 


23 


.07 


23 


.48 


1 


.48 


23 


25 


22 


.76 


23 


.47 


23. 


35 


1 


.03 


22 


.24 


21 


.93 


21.40 


S31 


22 


.94 


1 


.47 


23 


18 


22 


.82 


22 


.81 


23 


34 


1 


.37 


23 


.21 


22 


.64 


23 


.25 


23. 


35 


1 


.00 


22 


.33 


21 


.80 


21.61 


S32 


22 


.92 


1 


.49 


23 


24 


22 


.87 


22 


42 


23 


42 


1 


.32 


23 


11 


22 


65 


22 


.85 


23 


35 


1 


.03 


22 


.27 


21 


95 


21.58 



a The apparent magnitude is M — Z p — 2.51og (DN/sec). 

^The 5<t detection limit for a circular aperture with diameter of 3 times FWHM 
c Thc magnitude corresponding to the completeness 50% 
d The magnitude corresponding to the reliability 99% 



Table 4. B— , R— , and J— band merged catalog. 



ID 




RA 




Dec 






Bauto 


Error 


Barer 


Error 


Rauto 


Error 


Raper 


Error 


I AUTO 


Error 


Iaper 


Error 


Stellarity 


A 


B 


A 


n 


A 


1 


flags 


B01.00001 


18 


11 


06 


3 


+66 


■:>7 


11 


5 


22 


88 


0.12 


2.3 


00 


0.10 


99 


00 


99.00 


99 


00 


99.00 


99 


00 


99 


00 


99 


00 


99 


00 





89 





18 





11 





08 




B01.00002 


18 


11 


08 


9 


+66 


51 


35 


9 


21 


07 


0.10 


22 


35 


0.08 


99 


00 


99.00 


99 


00 


99.00 


99 


00 


99 


00 


99 


00 


99 


00 





00 





18 





11 





08 




EO 1.00003 


18 


11 


Oil 


3 


+67 


07 


40 


4 


23 


07 


0.16 


23 


32 


0.14 


99 


00 


99.00 


99 


00 


99.00 


99 


00 


99 


00 


99 


00 


99 


00 





03 





19 





12 





08 


ncar_bright_obj 


E01.00004 


18 


11 


09 


3 


+66 


51 


19 


9 


23 


59 


0.33 


24 


10 


0.27 


99 


00 


99.00 


99 


00 


99.00 


99 


00 


99 


00 


99 


00 


99 


00 





18 





18 





11 





08 


galaxy 


Trm aaaak 


18 


11 


09 


6 


+67 


07 


■:)4 


8 


23 


39 


0.25 


23 


96 




99 


00 


99.00 


99 


00 


99.00 


99 


00 


99 


00 


99 


00 


99 


00 





04 





19 










08 


ncar_brigfit_ob j 


E01.00006 


18 


11 


Oil 


6 


+66 


51 


22 


9 


23 


03 


0.21 


24 


02 


0.25 


99 


00 


99.00 


99 


00 


99.00 


99 


00 


99 


00 


99 


00 


99 


00 





09 





18 





11 





08 


galaxy 


E01.00007 


18 


11 


10 


2 


+67 


07 


32 





20 


85 


0.08 


22 


52 


0.09 


99 


00 


99.00 


99 


00 


99.00 


99 


00 


99 


00 


99 


00 


99 


00 





00 





19 





12 





08 


ncar_bright_obj 


E01.00008 


18 


11 


10 


5 


+67 


07 


51 


7 


21 


90 


0.11 


23 


19 


0.13 


99 


00 


99.00 


99 


00 


99.00 


99 


00 


99 


00 


99 


00 


99 


00 





01 





19 





12 





08 


ncar_bright_obj 


E01.00009 


18 


11 


11 


2 


+67 


05 


29 


6 


22 


82 


0.18 


23 


59 


0.18 


99 


00 


99.00 


99 


00 


99.00 


99 


00 


99 


00 


99 


00 


99 


00 





32 





19 





12 





08 




E01.00010 


18 


11 


11 


2 


+66 


50 


56 


2 


99 


00 


99.00 


99 


00 


99.00 


22 


15 


0.12 


22 


21 


0.12 


99 


00 


99 


00 


99 


00 


99 


00 





36 





18 





11 





08 


ncar_cdgc 


E01.00011 


18 


11 


11 


3 


+66 


50 


53 


4 


99 


00 


99.00 


99 


00 


99.00 


21 


91 


0.11 


22 


04 


0.12 


99 


00 


99 


00 


99 


00 


99 


00 





47 





18 





11 





08 


ncar_cdgc 


E01.00012 


18 


11 


11 


4 


+67 


05 


27 


2 


22 


35 


0.14 


23 


33 


0.15 


99 


00 


99.00 


99 


00 


99.00 


99 


00 


99 


00 


99 


00 


99 


00 





21 





19 





12 





08 


galaxy 


E01.00013 


18 


11 


11 


4 


+67 


07 


10 


7 


20 


55 


0.06 


20 


60 


0.06 


99 


00 


99.00 


99 


00 


99.00 


99 


00 


99 


00 


99 


00 


99 


00 





98 





19 





12 





08 


ncar_bright_obj 


E01.00014 


18 


11 


11 


4 


+66 


50 


51 


2 


99 


00 


99.00 


99 


00 


99.00 


22 


77 


0.16 


22 


27 


0.14 


99 


00 


99 


00 


99 


00 


99 


00 





46 





18 





11 





08 


ncar_cdgc 


E01.00015 


18 


11 


11 


5 


+66 


51 


19 


2 


99 


00 


99.00 


99 


00 


99.00 


22 


69 


0.17 


22 


46 


0.16 


99 


00 


99 


00 


99 


00 


99 


00 





43 





18 





11 





08 


ncar_cdgc 


E01.00016 


18 


11 


11 


6 


+67 


(15 


27 


1 


22 


95 


0.17 


23 


66 


0.19 


99 


00 


99.00 


99 


00 


99.00 


99 


00 


99 


00 


99 


00 


99 


00 





71 





19 





12 





08 




E01.00017 


18 


11 


11 


7 


+67 


07 


12 


8 


20 


50 


0.06 


20 


62 


0.06 


99 


00 


99.00 


99 


00 


99.00 


99 


00 


99 


00 


99 


00 


99 


00 





96 





19 





12 





08 


ncar_bright_obj 


E01.00018 


18 


11 


11 


7 


+66 


51 


04 


3 


99 


00 


99.00 


99 


00 


99.00 


22 


81 


0.19 


22 


53 


0.18 


99 


00 


99 


00 


99 


00 


99 


00 





40 





18 





11 





08 


ncar_cdgc 


E01.00019 


18 


11 


11 


I) 


+67 


07 


22 


9 


23 


59 


0.79 


23 


24 


0.14 


99 


00 


99.00 


99 


00 


99.00 


99 


00 


99 


00 


99 


00 


99 


00 





00 





19 





12 





08 


near_bright_obj 


E01.00020 


18 


11 


12 


1 


+66 


53 


23 


3 


99 


00 


99.00 


99 


00 


99.00 


18 


10 


0.05 


19 


19 


0.05 


99 


00 


99 


00 


99 


00 


99 


00 





03 





18 





11 





08 





CO 

oo 

Note. — (1) The identification number of the source. Each ID consists of the subfield name (Figure^ and number in the order of their RA. 
(2) and (3) The J2000.0 right ascension (RA) and the declination (Dec) of a source in a sexigesimal. 

(4)— (15) are the total magnitudes, its magnitude errors, aperture magnitudes and its magnitude errors of a source in the B—,R—, and J — band. The aperture magnitudes arc 

obtained with the aperture diameter of 3 times FWHM. The magnitudes are not corrected for the extinctions. 

(16) The stellarity index of a source from SExtractor. Sec Figurc[6]for the distribution of this index. 

(17—19) The galactic extinction values of B — , R— , and / — band from extinction map of Schlcgcl ct al. (1998). 

(20) The flags parameter. 'near_bright_obj ' for the spurious detection near bright stars, l ncar_cdgc' for the spurious detection near the edge of a single image, or 'galaxy' for the 
non-stellar sources classified at section |4.2| 



Table 5. Galaxy and stellar number counts. 



B — band R — band J — band 



mag 


galaxy 
number 


error 


corrected 
number 


star 
number 


galaxy 
number 




error 


corrected 
number 


star 
number 


galaxy 
number 




error 


corrected 
number 


star 
number 


15.75 


3 


.68B+00 


9.85E-01 


3.68E+00 


1 


46E+02 


1 


.13E+01 


1 


.73E+00 


1 


13E+01 


2 


.24E+02 


9 


21E+00 


1 


56E+00 


9 


21E+00 


2 


.72E+02 


16.25 


3 


■ 42B+00 


9.49E-01 


3.42E+00 


1 


74E+02 


1 


.71E+01 


2 


■ 12E+00 


1 


71E+01 


2 


.96E+02 


1 


.76E+01 


2 


15E+00 


1 


76E+01 


3 


■ 35E+02 


16.75 


6 


■ 58B+00 


1.32E+00 


6.58E+00 


2 


19E+02 


2 


■ 58E+01 


2 


61E+00 


2 


.58E+01 


3 


■ 27E+02 


3 


.30E+01 


4 


28E+00 


3 


31E+01 


4 


03E+02 


17.25 


1 


■ 26B+01 


1.82E+00 


1.26E+01 


2 


49E+02 


4 


.24E+01 


3 


34E+00 


4 


.24E+01 


3 


.91E+02 


5 


.93E+01 


5 


91E+00 


5 


95E+01 


5 


.04E+02 


17.75 


2 


■ 13B+01 


2.37E+00 


2.13E+01 


2 


92E+02 


8 


■ 58E+01 


4 


.75E+00 


8 


.58E+01 


4 


.58E+02 


1 


.03E+02 


9 


62E+00 


1 


03E+02 


5 


.84E+02 


18.25 


2 


97B+01 


2.80E+00 


2.97E+01 


3 


42E+02 


1 


■ 15E+02 


8 


.29E+00 


1 


■ 15E+02 


5 


■ 25E+02 


1 


.95E+02 


1 


42E+01 


1 


97E+02 


7 


06E+02 


18.75 


5 


.55B+01 


3.82E+00 


5.55E+01 


3 


73E+02 


2 


■ 19E+02 


1 


■ 27E+01 


2 


■ 20E+02 


5 


.99E+02 


3 


.34E+02 


3 


.08E+01 


3 


38E+02 


7 


.96E+02 


19.25 


8 


05B+01 


4.60E+00 


8.05E+01 


3 


91E+02 


3 


■ 73E+02 


2 


.81E+01 


3 


.78E+02 


7 


.24E+02 


5 


.07E+02 


7 


27E+01 


5 


30E+02 


8 


.81E+02 


19.75 


1 


■ 27B+02 


5.77E+00 


1.27E+02 


4 


78E+02 


6 


■ 00E+02 


5 


.98E+01 


6 


■ 20E+02 


8 


■ 28E+02 


7 


■ 29E+02 


1 


.02E+02 


7 


87E+02 


9 


.12E+02 


20.25 


1 


■ 99B+02 


2.62E+01 


2.02E+02 




38E+02 


9 


.25E+02 


9 


■ 37E+01 


9 


■ 78E+02 


9 


■ 36E+02 


1 


.08E+03 


1 


.24E+02 


1 


24E+03 


9 


■ 83E+02 


20.75 


3 


■ 76E+02 


5.74E+01 


3.89E+02 


6 


09E+02 


1 


■ 26E+03 


1 


■ 20E+02 


1 


■ 38E+03 


1 


.02E+03 


1 


■ 54E+03 


1 


49E+02 


2 


01E+03 


1 


■ 05E+03 


21.25 


6 


■ 59B+02 


8.96E+01 


7.10E+02 


7 


69E+02 


1 


■ 48E+03 


1 


■ 38E+02 


1 


.73E+03 


1 


00E+03 


2 


02E+03 


1 


70E+02 


3 


12E+03 


9 


■ 80E+02 


21.75 


9 


.65B+02 


1.12E+02 


1.11E+03 


8 


38E+02 


1 


.71E+03 


1 


■ 47E+02 


2 


21E+03 


8 


.55E+02 


















22.25 


1 


■ 41B+03 


1.35E+02 


1.88E+03 


9 


24E+02 


2 


36E+03 


1 


.71E+02 


3 


.37E+03 


7 


.24E+02 


















22.75 


2 


01B+03 


1.58E+02 


2.97E+03 


9 


31E+02 


3 


.37E+03 


2 


.02E+02 


5 


.56E+03 


6 


31E+02 


















23.25 


3 


.13B+03 


1.95E+02 


5.09E+03 


7 


14E+02 


6 


■ 12E+03 


2 


.73E+02 


1 


.12E+04 


3 


.98E+02 



















Note. — Units of the galaxy and stellar number densities and their errors arc Numbers 0.5mag 1 deg 



